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Beschreibung 

[0001] Die Erfindung bezieht sich auf ein Verfahren zur insgesamt treibstoffminimalen, rechnergestOtzten Ansteue- 
rung einer beiiebig festgelegten Anzahl und beliebig festgelegten Anordnung von DQsen an einem Raumfahrzeug. Die 

5 Dusen konnen dabei diskret oder auch kontinuierlich ansteuerbar sein. 

[0002] Fur die Ausubung von auBeren Kraften und Momenten an einem , Raumfahrzeug ist der Einsatz von Dusen 
in vieien SItuationen notwendlg. Da der Einsatz der Dusen Treibstoff verbraucht und dieser unmittelbar Elnflu3 auf das 
Startgewicht und die Lebensdauer von Raumfahrzeugen und insbesondere von Satelliten hat, kommt der Dusenan- 
steuerung mrt in der Summe minimalem Trelbstoffverbrauch eine wichtige Bedeutung zu. Diisenanordnungen, die 

10 prinzipieli einen niedrigen Treibstoffverbrauch emnogllchen, eine hohere oder gleiche Redundanz bei klelnerer DO- 
senanzahl aufweisen und/oder durch die Nutziast am Raumfahrzeug beeinfluBt sind, sind oft unregelm^Bige und 
schiefwinklige DQsenanordnungen. Aus diesen GrOnden ist die Forderung nach dertrelbstoffmininnalen DQsenansteue- 
rung, insbesondere auch fur beliebige (schiefwinkilge) feste DQsenanordnungen an Raunnfahrzeugen von besonderer 
Relevanz. 

15 [0003] Es ist bekannt, daB sich das Problem der treibstoffminimalen Dusenansteuerung auf ein iineares Optimie- 
rungsproblem zuniickfuhren Iaf5t, fur das der Simplexalgorithmus ein geeignetes und effizientes Losungsverfahren 
darstelit (siehe zum Beispiel RJ. Wiktor "IN^inimum Control Authority Plot: A Tooi for Designing Thruster Systems" In: 
Journal of Guidance, Control, and Dynamics, Vol. 17, No. 5, Sept.-Oct. 1994, pp. 4-5). Der Einsatz des Simplexalgo- 
rithmus scheitert jedoch haufjg an der benotigten Rechenzeit. So konnen oft die Echtzeitanforderungen zum Beispiel 

20 fur den Einsatz des Simplexalgorithmus Innerhalb des Regeltaktes des Lage- und Positionsregelkreises eines Raum- 
fahrzeuges nk:ht erfQIlt warden, weil die besonders in hohen Ori?its einsetzbaren Rechner den Simplexalgorithmus 
nicht ausrelchend schnell abariselten kdnnen. 

[0004] Eine Vielzahl der sich heute im Einsatz befindenden Raumfahrzeuge (Satelliten) besitzen eine Dusenanord- 
nung, bei der jede Duse zum kdrperfesten Koordinatensystem des Raumfahrzeugs zu einer Achse parallel und damit 

25 zu den zwei anderen Achsen orthogonal angeordnet ist. Diese orthogonale und achsenparallele Anordnung hat zur 
Folge, daB jede Duse nur Krafte in einem RIchtungssinn einer Achse nrichtung und Momente nur in einem Drehsinn 
um die beiden anderen Achsen erzeugen kann. Der Vorteii dieser Anordnung ist, daS die Steuerung elnes jeden Frei- 
heltsgrades des Raumfahrzeuges als stan'em Korper durch den Einsatz von jeweils nur zwei Dusen moglich Ist, ohne 
andere Freiheitsgrade zu beeinflussen. Die Ansteuerung der Dusen ist fur diese Dusenanordnung sehr einfach. Der 

30 Nachteil dieser DQsenanordnung besteht darin, daB mindestens zwolf DQsen notwendlg sind, um aile Freiheitsgrade 
gleichzeitig und beliebig steuem zu konnen. MuB der Ausfall von zwei beliebigen Dusen tolerierbar sein, sind 36 Dusen 
erforderlich. 

[0005] Demgegenuber liegt der Vorteii einer geelgneten schiefwinkligen, also nIcht orthogonalen und achsenparal- 
lelen DQsenanordnung im wesentlichen in folgenden drei Punkten: 

35 

1 . Bei frei wahibarer fester Dusenanordnung sind weniger als 12 Dusen erforderlich. 

2. Fur die Tolerierung des Ausfalls von zwei beliebigen Dusen sind erheblich weniger als 36 Dusen erforderlich. 

3. Bei treibstotfminimaler Ansteuerung der Dusen Ist ein gerlngerer Trelbstoffverbrauch moglich. 

40 [0008] Der Nachteil dieser DQsenanordnung besteht darin, daB die treibstoffmlnlnnale DQsenansteuerung die Losung 
eines linearen Optimierungsproblems erfordert und der Rechenautwandfur die Ldsung dieses Problems mit dem Sim- 
plexalgorithmus oft zu hoch ist. Aus diesem Grunde sind andere Verfahren entwickelt worden. die weniger Rechenzeit 
benotigen und durch die eine mehr oder weniger treibstoffminimleite, aber keine Insgesamt trelbstoffmlnimaie Dusen- 
ansteuerung erreicht werden kann. 

45 [0007] Beim sogenannten "Linear Programming" wird das DQsenansteuerungsproblem als Iineares Optimierungs- 
problem aufgefaBt und mit einschiSgigen Verfahren der linearen Optimierung, wie zum Beispiel den Simplexalgorith- 
mus, geldst. Der schon erwahnte Nachteil, der den Einsatz dieses Verfahrens in vtelen Fallen unm5giich macht, ist 
derzu hohe Rechenaufwand. 

[0008] Beim 'Table Look-up"-Verfahren werden alle treibstoffminimalen DQsenanordnungen fur EInheitskrafte und 
so -momente bezQglich einer Achse des kdrperfesten Koordinatensystenns des Raumfahrzeugs berechnet (Losen des 
linearen Optimierungsproblems) und in einerTabelle abgelegt. Sind im Einsatz Kr&fte und Momente bezQgltoh mehrerer 
Achsen gleichzeitig auszuQben, werden die entsprechenden Dusenansteuerungen bezuglich Kraft oder Moment und 
fur jede geforderte Achse und dem geforderten Richtungs-oder Drehsinn aus der Tabelle ausgelesen und mit dem 
entsprechenden Betragswert der Kraft- bzw. Momentenkomponente in dieser Richtung multipliziert. Den Ansteuerwert 
55 fur jede Duse erhalt man durch die Addition der einzelnen Werte fur die jeweilige Duse. Dieses Verfahren ist sehr 
einfach und sehr schnell, fuhrt im allgemeinen aberzu nicht treibstoffminimalen DQsenansteuerungen. 
[0009] Bei verfeinerten und erweiterten "Table Look-up" -Verfahren handelt es sich um Verfahren, die eine Look-up- 
Tabelle aniegen, die EintragefOrzusStzlich eingefuhrte Achsen enthalten. Durch diese Verfahren kdnnen treibstoffspa- 
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rende DOsenansteuemngen ennitteit werden. Zuglelch nimmt aber die Oberslchtlichkeit des Verfahrens ab und wegen 
derzunehmenden Gr6Be derTabellen der Rechenaufwand erheblich zu. 

[001 0] Andere Verfahren beruhen auf der Auswertung und Verbesseaing des Ergebnlsses, das mit dem "Table Look- 
up"-Verfahren ermittelt wird. Das Ziel dieser Auswertung besteht darin, besonders trelbstoffverschwendende DQsen- 
5 ansteuerungen zu verhindern bzw. geeignetzu korrigieren. Hier soli exemplarisch ein Verfahren, das auf einenn Fuzzy- 
Ansatz beruht, genannt werden (vgl. T.SuzuW, K.Yasuda, S.Yoshikawa, K.Yamada, N.Yoshida "An Application of Fuzzy 
Algorithm to Thruster Control Systems of Spacecraft" in: ISTS 94-C-1 6, 1 9th International Symposium on Space Tech- 
nology and Science, Yokohama. Japan, May 15-24, 1994, pp. 1 - 6). 

[001 1] Eine Treibstoffminimierung der DQsenansteuerung Insgesamt hat zum Ziel eine sok;he Ansteuerung der DQ- 
10 sen zu ermittein, die 

1. die geforderten Krafte und Momente am Raumfahrzeug ausubt, 

2. fur jede Duse einen nichtnegatlven AnsteuenYert iiefert, und 

3. die Summe aller Ansteuerwerte der Dusen minimiert. 

15 

[0012] Dies stellt - wie erwahnt - ein lineares Optimierungsproblem dar. Da das Losen des Problems in der Regel 
zu zeitaufwendig 1st, kann die Losung des quadratischen Optlmierungsproblems als Ersatzproblem durchgefuhrt war- 
den. 

[0013] Die Losung des Ersatzproblems fuhrt zur Ermittlung derjenigen DQsenansteuerung, die 

20 

1 . die geforderten Krafte und Momente am Raumfahrzeug ausubt, 

2. fur jede Duse einen nichtnegatlven AnsteuenYert Iiefert, und 

3. die Summe aller quadratischen Ansteuerwerte der Dusen minimiert. 

25 [0014] Dieses quadratische Optimierungsproblem kann uber die Bereitsteiiung der Moore- Penrose-invertierten Sy- 
stemmatrix der DQsenanordnung am Raumfahrzeug und eines Offsetvektors, der mit Hilfe der singularen Wertezerle- 
gung berechnet wird (vgl. z.B. B.Noble, I. W. Daniel "Applied linear algebra" Prentice-Hall International, INC. (1988), 
pp. 338 - 350), gelost werden. Die Losung des quadratischen Optlmierungsproblems kommt dann als Produkt der 
Moor-Penrose-invertierten Systemmatrix mit dem kommandlerten KrafterMomenten-Vektor und der Addition des Off- 

30 setvektors zustande. Die Addition des Offsetvektors bewlrkl, da(3 Im LSsungsvektor nur nicht negative Ansteuerwerte 
auftreten. Die Berechnung des Offsetvektors wird jedoch dann numerisch immer aufwendiger, wenn die Zahl der an- 
zusteuemden Dilsen urn zwei oder mehr uber der Dimension des kommandierten KrSfte-Momenten-Vektors liegL Im 
allgemeinen 1st die Losung nicht treibstoffHninimal, fallt aber mit der Losung des llnearen Optimierungs-problems zu- 
sammen, wenn die Zahl der anzusteuemden Dusen nur um eins groBer ist als die Dimension des kommandlerten 

35 Krafte-Momenten-Vektors (vergleiche. z.B. DE 195 19 732 A1) . Weiterhin Ist nachteilig, dal3 iangere Einschaltzelten 
Qberproportionai "bestrafT werden und daher seltener in den Ansteuerungen vorkommen. Fur den Betrieb der Dusen 
und ihrem Wirkungsgrad ist es jedoch gunstiger, sie mit langeren Einschaltzelten zu betreiben. 
[0015] Der Nachteil bestehender Verfahren besteht mithin darin, daB sie entweder nicht das lineare Optimierungs- 
problem losen und damit im allgemeinen nicht die treibstoffminimale Ansteuerung der Dusen bereitstellen oder aber, 

40 wenn bei den Ansteuerverfahren das lineare Optimierungsproblem gelost wird und der Simplexalgorithmus ausgefOhrt 
wird, dieses ubiicherweise zu rechenaufwendig ist. 

[001 6] Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der eingangs angegebenen Art bereitzustellen, mit 
dem das lineare Optimierungsproblem fur die treibstoffminimale Ansteuerung beliebig fester Dusenanordnungen zur 
gleichzeitigen Realisiemng von Kraften und Momenten an einem Raumfahrzeug in signifikant kurzerer Zeit gelost 

45 werden kann, als es bisher moglich isL 

[0017] Diese Aufgabe wird gemaS der Erfindung durch die im Anspruch 1 gekennzeichneten Merkmale gelost. 
[0018] Das Verfahren nach der Erfindung basiert also vorteilhaft auf einem. dem Simplexalgorithmus angelehnten 
Verfahren, das Jedoch aufgrund der konsequenten Ausnutzung der Besonderheiten be] der DQsenansteuerung bei 
Raumfahrzeugen und der Aufgabe der Verwendung des Qbllchen SImplextableaus zu elner RechenzeltverkQrzung von 

50 ca. einer GroBenordnung fuhrt. 

[001 9] Vortellhafte Ausgestaltungen und WefterfOhrungen des Verfahrens nach der Erfindung sind in den dem An- 
spruch 1 nachfolgenden Patentanspruchen gekennzeichnet. 
[0020] Das Verfahren nach der Erfindung basiert auf folgenden Grundlagen: 

S5 a) Das bei der Dusenansteuerung zu losende lineare Optimierungsproblem kann wie folgt aufgeschrieben werden: 

Gesucht ist der Ansteuervektor Kopi, der die Ansteuerwerte fOr alle DOsen enthSIt und folgende Bedingungen er- 
fiilien muB: 
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1. Durch J^pt mussen die geforderten KrSfte und Momente b^JJ^ am Raumfahrzeug ausgeubt werden. Das 
herBt, die Gleichung A^ys • j<opt ' emiWt werden, wobei die Systemmatrix ist, in der die drei Kraft- 
und die drei Momentenlcomponenten fur jede einzeine Duse in einer Spatte dieser Matrix eingetragen sind. 
Die Systemmatrix hat 6 Zeilen wegen der 6 mdglichen Freiheitsgrade in der Bewegung des Raumfahrzeugs, 
namlich 3 f<rafte in drei zueinanderorthogonalen Achsen und entsprechend 3 Momente urn diese Achsen und 
bei D Dusen D Spalten. 

2- kopt = [kopt 1 •• Kspti " KiptoF darf nur nichtnegative Werte enthaiten, 
d.h.: 



10 



ICop,,^0.j=1(1)D. 
3. die Summe alter Ansteuerwerte muB minimal sein (Zielfunktlon), d.h.: 



15 
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opt j 



= Min ! 
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b) Ein Simplextableau besteht aus einem Spaltenvektor der Basisvariablenindizes X^,, einem Zeilenvelctor der 
Nichtbasisvariablenindizes X^b, der Koefflzientenmatrix A^gb, dem Zeilenvelctor der auf die Nichtbasisvarlablen 
umgerechneten Zielfunktionskoeffizienten und dem Spaltenvelctbr der Basisvariablenwerte btat,: 
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Im Verfahren nach der Erfindung kommt das Simplextableau nicht In dieser vollen Fonn zum Einsatz. 

c) Fur die Umsetzung eines bellebigen 6-dimensionalen Krafle-Momenten-Vektors sind maximal 6 Dusen einer 
45 vorhandenen festen Dusenanordnung erforderlich Oder der kommandierte Krfifte-Momenten-Vektor ist mit der 

vorhandenen Dusenanordnung nicht reallsierbar. 

d) Fur jede Diisenauswahl I von 6 aus D Dusen existiertein BasisvariablenindizesvektorXb (I). Die Nummem oder 

Indizes (1 ,2 D) der in der Auswahl I enthattenen Dusen werden irn Basisvariablenlndizesvektor (i) zusam- 

mengefaSt und die Ansteuerwerte dieser Dusen sind die Werte der Basisvariablen des Simplextabieaus. Die Werte 

so der Basisvariablen sind entsprechend der Reihenfolge der Dusenindizes in (I) im Basisvariablenwertevektor 

btab(^ b'^om) abgelegt. Zujedem Basisvariablenlndizesvektor (I) gehort eine Basismatrix Bt, (I). Die Basismatrix 
wird gebildet, indem man diejenigen Spatten der Systemmatrix Agys, die zu den ausgewahlten Dusen gehoren, zu 
einer neuen Matrix spaltenweise zusammenstellt. Die Reihenfolge der Spalten in der Basismatrix entspricht der 
Reihenfolge der Indizes der Dusen im dazugehorigen Basisvariablenlndizesvektor Die Nichtbaslsmatrix B^b 

55 (1) wIrd aus denjenigen Spalten der Systemmatrix zusammengestellt, die zu den nicht in der Auswahl I enthaltenen 

Dusen gehoren. Ihre Indizes sind im Nichtbaslsvariablenindizesvektor ^ nb(l) zusammengefaBt. Die Reihenfolge 
der Spalten in der Nichtbaslsmatrix entspricht der Reihenfolge der Indizes der Dusen im dazugehorigen Nichtba- 
sisvariablenindizesvektor. Die Ansteuerwerte dieser DQsen sind Null. 
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e) 1st fur einen kommandlerten Krfifte-Momenten-Vektor bnom das Optlmaltableau Xfc^opt • ^.opt • Atab.opi. Etab.opt 
und btab opt gefunden, dann sind die Ansteuerwerte der DOsen, die durch die Indizes in Xnt opt adresslert werden, 
gleich Null und diejenigen Dusen, die den Basisvarlablen entsprechen, also durch die Indizes in adressiert 
werden, enthalten als Ansteuenwert den Wert aus der entsprechenden Zeile des optimalen Basisvariablenwerte- 

5 vektors bjab.opt- Basisvariablenindizesvektor Aj, gibt die optinnale Dusenauswahl an. 

f) Furden Fall der glelchzeltlgen Ausiibung von alien drel Kraftkomponenten und alien drel Momentenkomponenten 
1st der Einsatz von maximai 6 dervorhandenen DQsen erforderllch. Es gibt genau 

10 ^<^6 ) = 61 (D-6)l 

Moglichkelten (Kombinatlonen ohne Wiederholung), 6 aus D Dusen auszuwahlen. Werden die Indizes in den bei- 
den Vektoren der Basisvariablenindizes und Nichtbasisvariablenindizes, X^, und X^^, stets aufsteigend sortiert, so 
gibt es jeweils genau A(C^) dieser Vektoren. Fur eine teste Dusenanordnung gibt es (nur) eine bestimmte Anzahl 
IS von Kombinationen von Dusen, die die Dusenauswahl fur eine Losung des llnearen Optimlerungs-problems sein 

konnen. Diese Anzahl ist im allgemelnen deutlich kleiner als die Anzahl A(C|). 

[0021] Eine wesentliche Grundlage des Verfahrens besteht darin. daB stets mit sortlerten Basis- und Nichtbasisva- 
riablenlndizesvektoren gearbeitet wird. Im weiteren sind die Indizes oder Nummem der Dusen in den Basis- und Nicht- 
20 basisvariablenindizesvektoren z.B. stets der GroBe nach sortiert. Auf diese Weise wird es moglich, daB einer DQsen- 
auswahi genau ein Indizesvektor, genau eine Koeffizientenmatrix At^b = [^tab kjl« Qsnau ein Koeffizientenvektor der auf 
die Nichtbasisvariablen umgerechneten Zielfunktion und genau eine Basismatrix B^ und Nichtbasismatrix Bpb 
zugeordnet werden. 

[0022] Eine weitere wesentliche Grundlage des Verfahrens besteht darin, daB mit einem reduzierten Tableau gear- 
25 beitet wird. Das reduzierte Tableau besteht lediglich aus der Dusenauswahiidentifikationsnummer I und dem Basisva- 
riablenwertevektor bjab- Fur das Verfahren und die Arbeit mit dem reduzierten Tableau ist neben der Bereitstellung des 
Feides der Koeffizientenmatrizen aller Optimaltableaus und der PIvotspaltenlndizesmatrIx auch die Bereitstellung der 
DQsenauswahlidentiflkatlonsnummernmatrix, der Zeilenaustauschindizesmatrix und der Baslsindizesmatrix ertorder- 
lich. 

30 [0023] Das Verfahren besteht aus zwei Bestandtellen. Der erste, der Initlallsierungstell, dient zum einen der Bereit- 
stellung des Feides der Koeffizientenmatrizen alter Optimaltableaus 

.[0024] Die DQsenauswahiidentifikationsnummer lauft von 1 bis M ( I = 1 (1 )M ), wobei M die Anzahl alter Dusenaus- 
40 wahlmoglichkeiten ist, die fur einen belieblgen Krafte-IWomenten-Vektor b^om und die vorliegende beliebige und un- 
veranderliche Dusenanordnung optimal sind. Eine Koeffizientenmatrix eines Optimaltableaus iiegt im Sinne dieses 
Verfahrens dann vor, wenn aile Elemente im zur Koeffizientenmatrix gehorenden Vektor der auf die Nichtbasisvariablen 
umgerechneten Zielfunktionskoeffizienten 

Etab = Ptab (i^) "'^^^ negativ sind. Fur die Elemente der Koeffizientenmatrix muB daher gelten: 



45 



50 



' -J, Sab ic 3 - P tab j ^ 0 ""^ 3 - 1(1)6. 



[0025] Zum anderen werden im Initialisierungsteil des Verfahrens folgende AdreBmatrizen: die Pivotspaltenindizes- 
matrix P= [pi^], die Dusenauswahlidentifikationsnummernmatrix U = [UiJ, die Zeilenaustauschindizesmatrix T = [tjKl 
und die Baslsindizesmatrix L = PixJ bereitgestellt, damit auf die im zwe'rten Tell des Verfahrens benotigten Daten be- 
55 senders schnell zugegriffen werden kann und bei der Durchfuhrung der Basisaustauschschritte die Arbeit mit dem 
reduzierten Tableau moglich wird. Durch die Dusenauswahlidentifikationsnummer I = 1 ,2, .... M wird jeder Dusenaus- 
wahl, die optimal sein kann, und damit auch jeder Koeffizientenmatrix eines Optimaltableaus eine Identifikationsnum- 
mer zugewtesen. 
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[0026] In der PIvotspaltentndizesmatrix P = [PikI wird fur jede Zeile von jeder in F abgelegten Koeffizlentenmatrix 
der Index der dazugehorigen PIvotspalte s e {0, 1 , 2, 3 (D - 6)} abgelegt. Fur eine Koeffizlentenmatrix Ajat, und eine 
gewahlte Pivotzeile z ist die Pivotspalte s diejenige, fur die gilt: 



10 [0027] D.h. die zu einer Zeile einer Koeffizientennnatrix gehorige Pivotspalte ist diejenige Spalte, in der das Matrix- 
element der Koeffizlentenmatrix in dieser Zeile negativ ist und der Quotient aus dem Element des Velctors der auf die 
Nichtbasisvarlablen umgerechneten Zlelfunktionskoeffizienten In der Spalte und dem Betrag des entsprechenden Ma- 
trixelements aus der Koeffizientenmatrix am kleinsten ist. 

[0028] IVIlt Hiife der Pivotspaftenindlzesmatrix Ist fur eine Dusenauswaiil I, fur die In F eindeutig die dazugehorige 
15 Koeffizientenmatrix A^{\) abgelegt Ist, und gewahlter Pivotzeile z die Angabe der entsprechenden Pivotspalte mit s 
(AtabC). z) = P iz der Pivotspaltenlndizesmatrix mogllch, ohne dal3 dazu die Koeffizientenmatrix A]^{\) Oder der 
Vektor der auf die Nichtbasisvarlablen umgerechneten Zlelfunktionskoeffizienten £tab(l) bendtigt wird. 
[0029] Existiert zu einer bestimmten Zeile keine Pivotspalte, so ist das in der Pivotspaltenlndizesmatrix durch eine 
Null (p iz = 0) angegeben. 

20 [0030] In der Dusenauswahlidentifikationsnummernmatrix U = [u,J istfur jede DQsenauswahlidentlfikationsnummer 
I = 1 ,2, .... M und jede Zeile z = 1 ,2,..., 6 jeder. in F abgelegten Koeffizientenmatrix die Dusenauswahlidentifikations- 
nummerderjenigen Dusenauswahl angegeben, die nach dem Basisaustauschschrittvoriiegt, wenn vordem Basisaus- 
tauschschritt die Dusenauswahl I voriag und die Zeile z als Pivotzeile fiir den Basisaustauschschritt gewahit wurde. 
D.h., wenn zum Zeitpunkt n die Dusenauswahl l[n] voriiegt und ein Basisaustauschschritt mit der Pivotzeile z[n] durch- 

25 gefOhrt wird, kann die Identlfikatlonsnummer der neu entstehenden DQsenauswahl ohne Verwendung von Atab(l) nnit 
l[n + 1] = u ijnjjjnj direkt aus der DQsenauswahlidentifikationsnummemmatrix ausgelesen werden. Dies wird dadurch 
m6glich, daB durch die Kenntnls des aktuellen Tableaus (bzw. seiner Identlfikatlonsnummer I) und der Wahl der Pivot- 
zeile die Pivotspalte impllzlt festgelegt Ist. Zur Ermlttlung des Folgetableaus Ist also die expllzite Kenntnls der Pivots- 
palte s[n] unnotig. 

30 [0031] In der Zeilenaustauschindizesmatrix T = [t,J ist fur jede Zeile z = 1 ,2, .... 6 von jeder der in F abgelegten 
Koefflzientenmatrizen der Zeilenindex des Pivotelements im (sortlerten) reduzierten Tableau angegeben, das nach. 
dem Basisaustauschschritt voriiegt, wenn vor dem Basisaustauschschritt die Dusenauswahl I voriag und die Zeile z 
als Pivotzeile gewahit wurde. 

[0032] In der Basisindizesmatrix L = Q,] ist fur jede Dusenauswahlidentifikationsnummer der Vektor der Indizes der 
35 in der Dusenauswahl i enthaltenen DQsen, also der zu I geh5rige Basisvariablenlndizesvektor X||,(l) = 1 angegeben. 
Die DQsenindlzes in den Vektoren l| sind der GroBe nach sortiert. 

[0033] Durch den Initiaiislerungsteil des Verfahrens wird weiterhin eine Start-Dusenauswahlidentifikatlonsnummer I 
[Oj = und die dazugeh6rlge Invertierte Baslsmatrix 

^ irt= (SbCstart))'"" bereltgestellt. Aus numerischen Grunden werden alle Werte In F und Bj^^, die einen sehr kleinen 

40 Betrag e, (z. B.: e = 1 0'®) haben, zu Null Null gesetzt. 

[0034] Im zweiten Tell des Verfahrens wird zur Realislerung der Im. Krafte-IVIomentenvektor kommandierten Krafte 
und Momente die optimale Dusenansteuerung, d.h. die Ansteuerwerte der einzuschaltenden Dusen ermittett. Nur die- 
ser Tell des Verfahrens muB den gestellten Echtzeitanforderungen genugen. Dieser zweite Tell besteht wiedemm aus 
zwei nacheinander auszufuhrenden Teilen. Zuerst wird ein dual zulassiges reduziertes Tableau als Starttableau enmlt- 

45 telt Im zweiten Tell wird die erforderiiche Anzahl von Basisaustauschschritten durchgefuhrt. 

[0035] Ein dual zulassiges reduziertes Staritabieau wird wie foigt emfilttelt: Fur die Bereitstellung eines dual zutas- 
sigen reduzierten Tableaus ist ledlglich die Transformation des aktuellen Krfifte-Momente-Vektors b^m auf die fur den 
Start der Basisaustauschschritte festgelegte DQsenauswahl l[0] = \^ erforderitoh. Dazu dient die schon bereitgesteilte 
invertierte Baslsmatrix B^^^, und es gilt: 



[0036] Mit der schon bereitgestellten Dusenauswahlidentifikationsnummer fur den Start l[0] = Igtart und dem neu 
55 berechneten Vektor bt^b [0] = steht ein reduziertes und dual zulassiges Starttableau [l[0] [0]] bereit. 

[0037] Dann erfolgt die AusfQhrung der Basisaustauschschritte (zyklischer Verfahrensteil): 
Der folgend beschriebene Verfahrensteil muB wiederholt solange ausgefuhrt werden, bis die optimale Dusenauswahl 
mit den optimalen Ansteuerwerten gefunden wurde Oder festgestellt wird, daB es fQr den aktuell voriiegenden KrSfte- 
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, mit j = 1{1){D - 6). 
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Momenten-Vektor und die vorliegende Dusenanordnung keine Ldsung gibt. 

[0038] Da von einem dual zulassigen Simplextableau In Fomn des reduzlerten Starttableaus ausgegangen wird und 
durch die Basisaustauschschritte immer wieder nur eine Dusenauswahl erzeugt werden kann, deren Simplextableau 
dual zulassig ist, kann be! der Durchfuhrung der Basisaustauschschritte immer wieder auf die im Initialisierungsteil 
5 des Verfahrens bereitgestellten Daten insbesondere im Feld der Koeff izientenmatrizen aller optimaien und damit gleich- 
zeitig aller dual zulassigen Simplextabieaus und in der Pivotspaltenindizesmatrix zumckgegriffen werden. 

Erster Schritt: 

10 [0039] Finde einen negativen Wert In b^[n]. n In eckigen Klammern bezelchnet die Nummer des Iterationsschrittes. 
Enthaitbjab[n] keinen negativen Wert, llegt die Losung des linearen Optimlerungsproblems vor, und es kann der Vektor 
der Dusenansteuerwerte kgpt autgebaut werden. Enthalt bt^b [n] negative Werte, dann wird eine dieserZeilen als Pivot- 
zelle z[n] gewahlt. Es ist giinstig, die Zeile des betragsgroBten negativen Elements In bta5[nl als Pivotzeile zu wahlen. 

15 Zweiter Schritt: 

[0040] Mit Kenntnis der Dusenauswahlidentifikationsnummer l[n], und der im ersten Schritt gewahlten Pivotzeile z 
[n] wird aus der Pivotspaltenindizesmatrix P der dazugehorige Pivotspaltenindex mit s[n] = p i^njzjnj ausgelesen. Ist In 
P kein Pivotspaltenindex angegeben (p i^n^^^] = 0). e'"® andere Pivotzeile ausgewShlt werden. Ist zu keiner der 
20 moglichen PIvotzellen In der Pivotspaltenindizesmatrix eine Pivotspalte Indizlert, existlert kelne LSsung des linearen 
Optimlerungsproblems. Sind Pivotzeile und Pivotspalte ennittelt, wirdzum dritten Schritt Qbergegangen. 
[0041] Im dritten und Im vierten Schritt wird der eigentliche Basisaustauschschrrtt durch Ennittlung des nSchsten 
reduzlerten Simplextabieaus ausgefuhrt: 

25 Dritter Schritt: 

[0042] Mit Hllfe der Koeffizlenten In der Pivotspalte s[n] von Atab[n] wird aus b^n] der neue Basisvariablenwerte- 
vektor btab neu berechnet. Die benotigten Elemente des Plvotspaltenvektors a aus fi^M konnen direkt aus dem Feld 
der Koefflzientenmatrizen aller Optimaltableaus mit Bt^b k8[niM = ^s[n]i[nj ausgelesen werden. 
30 Die Elemente des neuen Baslsvarlablenwertevektors werden wie folgt berechnet: .... 
Der neue Wert In der Pivotzeile z[n] wird berechnet nach: 

*^iat nou z[ni " f ' 

35 nn]s[n]l[n] 



40 



wobei f2[n]s[n] i[n] ^as entsprechend Indizierte Pivotelement aus Aj^b [n] ist. 
[0043] Die neuen Werte in den anderen Zellen werden berechnet nach: 



wobei: (/c = 1 (1 )6) aber {k ^ 2[n]) 
[0044] Da der neuberechnete Basisvariablenwertevektor bjab.neu ^uiji fiktiv angenommenen neuen sortierten Basis- 

45 variablenindizesvektor X^[r\ + 1] passen mu3, ist im allgemelnen der neuberechnete Baslsvariablenwertevektor ent- 
sprechend umzuordnen. Dazu sind im neuberechneten Baslsvariablenwertevektor btg^.nw ledlglich die Werte zwischen 
der aktuellen Pivotzeile z[n] und der mit Hllfe von l[n] und z[n] aus der Zeilenaustauschindizesmatrix T = [t, J auszule- 
senden Austauschzelle z^usch W = * i[n]2[ni zyklisch zu tauschen. Die Austauschzelle ZjauscJ"] 'st die Zeile, wohin der 
Wert des neuberechneten Basisvariablenwertevektors plaziert wird, der In btab ^eu Pivotzeile z[n] steht. Nach 

50 diesem Umordnen des neuberechneten Basisvariablenwertevektors llegt der fOr den f olgenden Basisaustauschschritt 
benotigte Baslsvariablenwertevektor bjabCn + 1] vor. 

Vlerter Schritt: 

55 [0045] l^it Kenntnis von l[n] und der Pivotzeile z[n] wird mit l[n + 1] = Uj^nj^jn] die Dusenauswahlldentlfikationsnummer 
ermitteit, diezum neuberechneten und nach aufsteigenden Basisvariablenindizes sortierten Baslsvariablenwertevektor 
bjgj, [n + 1] gehort. Der Basisaustauschschritt ist damit. vollzogen, das neue reduzierte Tableau [l[n + 1] bj^b [n + 1)1 
Hegt vor, und es kann die Ldsung des linearen Optlmierungs-problems mit dem ersten Schritt des zyklischen Verfah- 
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rensteils fortgesetzt werden. 

[0046] Bricht das Verfahren mit einem Basisvariablenwertevektor, der keine negativen Werte mehr enthalt. ab, so 
liegt der optimale 6-dimensionale Basisvariablenwertevektor btab.opt = [%ndi vor, und es kann mit HItfe von l[nendl 
BUS der Basisindizesmatrix L derzum aktuelien Basisvariablenwertevektor gehorende Vektor der Basisvariablenindizes 
5 HDit XbCRendl = i[nendl ausgeiesen werden. Mit Kenntnis der belden Vektoren bjaJnenJ ""^ ^"enJ ^ird der optimale 
D-dimenslonale Ansteuervektor fOr die Dusen k^p^ aufgebaut. 

[0047] FQr die Losung des Ilnearen Optimierungsprobiems fQr unterschledltehe Krafte-Momenten-Vektoren sind im 
allgemeinen auch unterschiedlicfi viele Basisaustauschschritte n^j^ e {0, 1, .... n^ax) erforderlicli. Bei Simulations- 
rechnungen mit einer Dusenanordnung von 12 Dusen wurden in sehrseltenen Fallen bis zu n^^ax^ ^ Basisaustauscli- 

10 schritte durchgefuhrt, obwohi im Durchschnitt nur ca. 3,5 Basisaustauschschritte notig waren. Bei geeigneter Wahl 
des Starttabieaus [Igtart ^tab.start] ^ind jedoch nur in ca. 0,5 % aller Falle mehr als 5 Austauschschritte notwendig. Im 
allgemeinen bleibt es ohne merkliche Konsequenzen fur die Anwendung, wenn die Anzahl der maximal auszufuhren- 
den Austauschschritte auf einen Wert begrenzt wird, bei dem nahezu alle Falle vollstandig bearbeitet werden und in 
den Fallen, In denen doch noch negative Werte Im Basisvariablenwertevektor stehen. diese zu Null gesetzt werden. 

IS Auf diese Welse kann der Worst-Case-Rechenaufwand welter gesenkt werden. Anwendungen dieser Verfahrensmdg- 
lichkeit bei der Simulation in geschlossenen Regelkrelsen zelgten keine bemerkbare Verschiechterung der Regelgute 
Oder eine Erhohung des Trelbstoffverbrauchs. 

[0048] Enthalt der Basisvariablenwertevektor nach der Ausfuhrung der maximal eriaubten Anzahl von Baslsau- 
tauschschritten noch negative Werte, istes prinzipiell auch moglich, daB keine Losung des Ilnearen Optimierungspro- 
20 blems existiert. Dieser Fall kann dann auftreten, wenn eIne Dusenanordnung vorliegt, die nicht bellebige Krafte und 
Momente aus dem Bereich aller zugelassenenen Krafte-Momenten-Vektoren erzeugen kann. Damit In diesen Fallen 
keine zu groBen Fehlkrafte und Fehlmomente erzeugt werden, kann in Abhanglgkeit vom Betrag der negtlven Elemente 
im Basisvariablenwertevektor entschleden werden, ob die negativen Elemente zu Null gesetzt werden und der An- 
steuervektor fur die Dusen aufgebaut wird oder generell keine Ansteuerung der Dusen erfolgt. 
25 [0049] Kommen bei der Abarbeltung des SImplexalgorfthmus entartete Optlmaltableaus vor, dann ist es theoretisch 
moglich, dal3 der Algorithmus in einen Endloszyklus lauft und nicht konvergiert. In Simulatlonsrechnungen mit unter- 
schiedllchen DQsenanordnungen Ist dieser Fall alierdings noch nie eingetreten und wird deshalb auch nicht explizit 
berucksichtigt Fur den dualen SImplexalgorithmus kann man durch die gewahlte Dusenanordnung sicherstellen, daB 
keine dual entarteten Optlmaltableaus existieren. Fur den primalen SImplexalgorithmus, also das Auftreten primal 
■ 30 • entarteter Optlmaltableaus, ist das so nicht moglich. 

[0050] Soil bei der Ermittiung der treibstoffminimalen Triebwerksansteuerung der unterschiediiche spezif ische Treib- 
stoffverisrauch C|, j=1 (1)D der etnzeinen Triebwerke berOckslchtlgt werden, liegt die Zielfunktlon 
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[0051 ] Die Berucksichtigung der Gewichtskoefflzienten Ist In dem Verfahren fQr den on-line auszufOhrenden Verfah- 
rensteii ohne Mehraufwand mdglch. 

[0052] Ffiiit eIne Duse aus, §ndert sich die Dusenanordnung im SInne dieses Verfahrens. Das Verfahren auf die 
neue Dusenanordnung anzuwenden erfordert, daB aile Daten, die im Initialisierungsteii des Verfahrens bereitgestellt 

45 wurden, fQr die neue DQsenanordnung bereitgestellt werden. 

[0053] Da sich die neue Dusenanordnung nur durch den Wegfall einer DQse von der alten DQsenanordnung unter- 
scheldet, gibt es Mogiichkeiten, die neuen Daten durch eine Aktualisierung der alten Daten zu emiitleln. Durch die 
Bereitsteilung eines erweiterten Feldes von Koeffizlentenmatrlzen, in dem auch diejenigen Koeffizientenmatrizen ab- 
gelegt sind. die durch den Ausfall von Dusen Koeffizientenmatrizen von Optlmaltableaus werden, und entsprechend 

50 erweiterter AdreBmatrizen, kann der Aktualislerungsaufwand stark reduzlert werden und fQr bestimmte Matrizen ganz 
entf alien. 

[0054] SchlieBlich besteht die Moglichkelt, die Daten fur alle DQsenanordnungen, die durch den Ausfall einer Duse 
zu erwarten sind, vorab bereitzustellen, und bei Ausfall einer bestimmten Duse dann auf die entsprechenden neuen 
Daten umzuschalten. Bis zum Ausfall einer nachsten DQse mussen fur die dann zu erwartenden DQsenanordnungen 
55 die Daten wiederum bereitgestellt werden, um dann durch Umschalten auf den nSchsten DQsenausfall reagieren zu 
konnen. 

[0055] Auch bei einer Verschiebung des Schwerpunktes des Raumfahrzeuges kommt es zu einer Veranderung der 
Dusenanordnung im Sinne dieses Verfahrens. Im allgemeinen f uhrt das zu einer VerSnderung der Systemmatrlx. Daher 
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mussen alle von ihr abgeleiteten Daten neu bereitgestellt werden. 

[0056] Das zuvor beschriebene Verfahren kann zur weiteren Einsparung von Rechenaufwand und Speicherplatz- 
bedarf noch modifiziert werden: 

5 - In voranstehender Beschreibung wurde stets von der Steuerung aller 6 Freiheitsgrade des Raumfahrzeugs aus- 
gegangen, Das Verfahren ist auch auf die Steuerung von 1 bis 5 Freiheits-, graden anwendbar. Es ist ledigiich die 
Streichung bzw. die ISIichtberQclcsichtigung der den nicht zu steuemden Freiheitsgraden entsprechenden Zeilen 
In der Systemmatrix und der von ihr abgeleiteten Velcloren und Matrizen erforderlich. Desgleichen ist die Steuerung 
einer sich st&ndig andemden Teiimenge aller 6 Freiheitsgrade moglich. Liegt eine schlefe DQsenanordnung vor, 

10 blelben Jedoch die nicht gesteuerten Freiheitsgrade im allgemeinen nicht unbeeinfluBt. 

IVIussen am Raumfahrzeug nicht beliebige Krafte und Momente ausgeubt werden, ist der Einsatz einer reduzierten 
DQsenanordnung moglich. Das vorllegende Verfahren Ist auch fur derartige Anordnungen anwendbar. Es sollte 
nur gewahrleistetsein, daB die kommandierten Krafte und Momente des Kommandos durch die vorllegende Du- 
15 senanordnung realisierbar sind. Ist das Kommando nicht realisierbar, gibt das Verfahren an, daB keine Losung 
des linearen Optimierungsproblems existiert. 

In Abhangigkeit von der vorllegenden DQsenanordnung kommt es vor, daB von den im Feld der Koeffizientenma- 
trizen aller Optlmaltableaus abgelegten Koeffizientenmatrizen bestimmte Spalten nie als PIvotspalte ausgewahlt 

20 werden. Nach Berechnung der Pivotspaltenlndlzesmatrix P konnen die nicht benotlgten Spalten aller Koeffizien- 

tenmatrizen einfach emilttelt werden. Es Ist moglich, ein Feld aller relvanten Pivotspaltenvektoren V zu erstellen, 
In dem nur die Spalten der Koeffizientenmatrizen aller Optimaltableaus abgelegt sind. die tatsachlich als Pivot- 
spaltenvektoren benotigt werden. Dadurch kann eine signifikante Spelcherplatzerspamis erreicht werden. Die 
Adressierung der Matrixspalten im Feld aller relevanten Pivotspaltenvektoren erfordert keinen erhohten Rechen- 

25 aufwand. Jede Spaite hat Im Feld aller reievan-ten Pivotspaltenvektoren ihren Spaltenindex. Wird In einer veran- 

derten Pivotspaltenindizesmatrix, der Pivotspalteniden-tifikationsnummemmatrix 0 = [OnJ nicht der Spaltenindex 
eines Plvotspaltenvektors In einer Koefflzlentenmatrix angegeben, sondem glelch die Pivotspaltenidentlfikations- 
nummer als Spaltenindex r im Feld aller relevanten Pivotspaltenvektoren V, dann kann Im zweiten Schritt des 
zyklischen Verfahrensteils mit r = o„j gleich der Spaltenindex Im Feld aller relevanten Pivot-spaltenvektoren V = 
[Vkrl angegeben werden. Die Enmittlung der Elemente des neuen Baslsvariablenwertevektors im dritten Schritt des 
zyklischen Verfahrensteils andert sich dann wie folgt: 

35 ^nn]iln] 

wobel V2[n] f^n] entsprechend indizierte Pivotelement aus Atab[n] 1st. Die neuen Werte in den anderen Zeilen 
werden berechnet nach: 

40 

wobei:{k= a&er(/r?&2[n]) Die weiteren Verfahrensbestandteile bleiben unverandert. 

45 - Llegt ein Feld aller relvanten Pivotspaltenvektoren V vor, kann Insbesondere bel Dusenanordnungen, dl6 bestimm- 
te Symmetrien aufwelsen, festgestellt werden, daB relevante Pivotspaltenvektoren unterschledlicher Optimalta- 
bleaus numerisch identisch sind. Da das mehrmallge Ablegen glelcher Pivotspaltenvektoren nicht notwendig ist, 
kann erneut weiterer Spelcherplatz eingespart werden. Beim sukzesslven Aufbau des Feldes aller relevanten 
Pivotspaltenvektoren im Initialisierungsteil des Verfahrens kann uberpruft werden, ob ein numerisch identischer 

so Pivotspaltenvektor zum gerade ennittelten Pivotspaltenvektor schon in V abgelegt wurde. Ist dies der Fall, kann 

das erneute Ablegen dieses relevanten Pivotspaltenvektors entfallen und In der zu V gehorigen Pivotspalteniden- 
tifikatlonsnummem matrix 0 wird die Pivotspaltenlndentifikationsnummer (Spaltenindex) des schon Im Feld aller 
relevanten Pivotspaltenvektoren abgelegten Pivotspaltenvektors eingetragen, 

55 - Besteht das Erfordernis, den Spelcherbedarf des Verfahrens noch welter zu reduzieren, kann eine Verfahrensva- 
riante eingesetzt werden, bel der es nicht erforderlich Ist, alle relevanten Ph^otspaltenvektoren, sondem nur noch 
Reprasen- , tantenvektoren im Rechnerspeicherabzulegen. Reprasentantenvektoren sind relevante Pivotspalten- 
vektoren, aus denen mit nur minlmalem Rechenaufwand weitere relevante Prvotspaltenvekloren berechnet werden 
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konnen. 

[0057] Jeder Pivotspaltenvektor a = [a^ag a] bl^t gehort zu einer bestimmten Basis, z.B. 
und zu eIner Nichtbasisvariablen, z.B. = V Wenn aj verschleden von Null ist, laSt sich bus 
derjenige Pivotspaltenvektor 

?neu ^ §('2. K) = fa„eu.1 ^neui ^neu.!,]^ 

mit 



a; 

berechnen, der zur Basis 

und zur Nichtbasisvariablen = A, gehSrl. So kann aus a(l^, Xr) n^^it der Angabe des Austauschlndexes i der neue 
Pivotspaltenvektor a(l2, A,) berechnet werden. . . 

[0058] Auf diese Welse kdnnen aus einem z.B. sechsdimenslonalen Pivotspaltenvektor bis zu sechs wettere Plvot- 
spaltenvektoren berechnet werden. Fur die Anwendung im Verfahren zurtrefbstoffinlnlnrialen DQsenansteuerung kann 
man sich jetzt auf das Speichem derjenigen relevanten Pivotspaltenvektoren als Reprasentantenvektoren beschran- 
ken, die erforderlich sind, unn alle relevanten Pivotspaltenvektoren wie angegeben zu berechnen. 
[0059] Fur die Emriittlung einer nnoglichst klelnen Anzahi von Reprasentantenvektoren lassen sich sehr gunstig Ge- 
netische Algorlthnnen als Optiomlerungsverfahren einsetzen. 

[0060] Fur eine geeignete Adresslerung der Reprasentantenvektoren in der Matrix aller Reprasentanten sind statt 
der Pivotspaltenindizes - bzw. Pivotspaltenindentifikatlonsnumnnemmatrixeinmal eine Reprasentantenldentifikations- 
nunnmernmatrlx R und welterhin eine AustauschindlzesmatrlxS erforderlich. In der Reprasentantenidentifikatlonsnum- 
nriemmatrix R wird fOr jede DQsenauswahlidentlflkationsnummer und jede Plvotzeile die Identlflkatlonsnummer des 
Reprasentantenvektors angegeben, derfOrdle Berechnung des tatsachiich bendtigten PIvotspaltenvektorsverwendet 
werden kann. In der AustauschlndizesmatrixS wird fur jede Dusenauswahlidentifikatlonsnumnner und Jede Plvotzeile 
der Zeilenlndex des Austauschelements angegeben, urn aus dem Reprasentantenvektor den tatsachiich benotigten 
Pivotspaltenvektor berechnen zu konnen. Stehtan einer Stelle in der Austauschindexmatrix eine Null, bedeutet dies, 
da3 der adressierte Reprasentantenvektor unmlttelbar der erforderliche Pivotspaltenvektor Ist und daher die Unrirech- 
nung in den aktuell benotigten Pivotspaltenvektor entfallt. 

[0061] Da stets mit nach Basisvariablenindrzes sortierten Pivotspaltenvektoren gearbeitet wird, muB der zunachst 
aus dem Reprasentantenvektor berechnete Pivotspaltenvektor noch umsortiert werden (z.B. nach aufsteigenden Ba- 
slsvariablenindizes), unn benutzt werden zu kdnnen. Wird die Sortierung derGroBeder Basisvariablen vorgenommen. 
gibt es bel eInem sechselementlgen Vektor genau 30 Mdgllchkeiten der Umsortierung, urn die erfor-derllche Ordnung 
wieder herzustellen, Bei einer erforderiichen Umsortierung wird stets ein Vektorelement an einer Position im Vektor 
entnommen und auf eine andere Position geschrieben. Die Vektorelemente zwischen den beiden Positionen werden 
entspre-chend verschoben. 

[0062] In einerSortieridentifikationsnummernmatrixG = [g , J wirdfurjeden Reprasentanten und jede Austauschzelle 
die Sortieridentifikationsnummer angegeben, nach der der Vektor umgeordnet werden muB, urn den benotigten Pivot- 
spaltenvektor zu erhalten. 

[0063] Aus einem Reprasentantenvektor konnen durch Baslsaustausch und/oder oben beschriebene Umordnung 
bis zu 7 x 31 = 21 7 Vektoren gebildet werden. Werden aite diese Mogiichkeiten berOcksichtigt, kann es aufgrund von 
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Symmetrien in derTriebwerksanordnung vorkommen, daBweniger ReprSsentantenvektoren erforderlich sind, ais wDr- 
den nurdiesieben aus einem Reprasentantenvektorgebildeten relevanten Pivotspattenvektoren In Betracht gezogen. 
[0064] Eine weitere Verfahrensvariante ergibt sich daraus, daB auf den Basisaustausch am Reprasentantenvektor 
verzichtetwirdund nurdieSO Umordnungsmoglichkeiten des Reprasentantenvektorsberucksichtigtwerden. Auf diese 
Weise reprasentiert ein Reprasentantenvektor 31 Vektoren. Aufgmnd von Symmetrien in derTriebwerksanordnung 
konnen auch auf dIese Weise mehrere relevante Pivotspaltenvektoren generiert werden. Im Vergleich zur letzten Ver- 
fahrensvariante warden im allgemeinen jedoch mehr ReprSsentantenvektoren benotigt, jedoch entfallt gleichzeitig der 
Rechenaufwand fur den Basisaustausch am Reprasentantenvektor, und es ist die Austauschindizesmatrix S nicht 
erforderlich. 



PatentansprOche 

1. Verfahren zur insgesamt treibstoffminimalen rechnergestutzten Ansteuerung einer beliebig festgelegten Anzahl 
und beliebig festgelegten Anordnung von Dusen an einem Raumfahrzeug, 
dadurch gekennzeichnet, 

daB ein treibstoffminlmaler Ansteuervektor (k^pi) mit ausschlieSlich nIcht negativen Ansteuerwerten fur alle 
vorgegebenen Dusen mit Hilfe des dualen Simplexalgorithmus gebildet wird, 

daB dazu jeder Dusenauswahl aus der Gesamtzahl von DQsen am Raumfahrzeug, die die optimale DDsen- 
auswahl fur die treibstoffminimale Realisierung eines beliebigen Kommandovektors i^om) darstellt, eine DQ- 
senauswahiidentiflkatlonsnummer (1) zugeordnet wird, 

da(3 fur die Emnittlung der optimalen Ansteuerung zunachst der Kommandovektor (bj^o^) der die auszuubende 
Kraft und/oderMomentenkomponenten enthaft, durch Multiplikation mit einer Transformatlonsmatrix {(Bsja^rt)'^) 
in ein en Startvektor (b=bstart) transfonnlert wird, 

daB danach der duale Simplexalgorithmus durchgefuhrt wird, bis der Vektor der Ansteuenverte (b) nur nicht- 
negative Ansteuerwerte enthStt oder aber die NIchtexistenz einer Losung angezeigt wird, 
und daB dabei der duale Simplexalgorithmus realisiert wird, Indem vorab alle Oaten, die fQr die Abarbeltung 
des dualen Simplexalgorithmus notwendig sind und sich off-line (vorab) berechnen lassen, emriittelt werden, 
in Datenfeldern Im Rechnerspelcher abgelegt werden und beim Etnsatz des dualen Slmplexatgorithmus w&h- 
rend des on-line-Betriebes durch geeignete Adressierungen aus dem Rechnerspelcher abgerufen und in der 
Rechnung verwendet werden, wobei nur mit einem reduzlerten SImpIextableau gearbeitet wird, indem lediglich 
die DQsenauswahlidentifikationsnummer und der dazugehdrige Vektor der Ansteuenverte (b) gefuhrt werden. 

2. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, 

daB die Transfomriationsmatrix ((B^))*"*) die invertierte Basismatrix derjenigen Dusen darsteiit, die zur Startdu- 
senauswahl geh&ren, diese Dusenauswahl fQr einen beliebigen Kommandovektor optimal ist, diese DQsenauswahi 
eine Start-Dusenauswahlidentifikationsnummer (I^tart} zugeordnet ist und der Startvektor (b==b^rt) Ansteuer- 
werte dieser Dusen enthalt. 

3. Verfahren nach einem der Anspruch e 1 oder 2, 

dadurch gekennzeichnet, 

daB jede Duse am Raumfahrzeug eine Identifikationsnummer erhalt, und die Werte im Vektor der Ansteuer- 
werte fur eine bestlmmte Dusenauswahl stets In einer festgelegten Ordnung anhand der Identifikatlonsnum- 
mem der Dusen sortlert werden, 

daB bei der Abarbeitung des dualen Simplexalgorithmus nach der Ermittlung einer Pivotzeilennummer und 
mit Kenntnis der DQsenauswahlidentiflkatlonsnummer der aktueil vorllegenden DQsenauswahi die entspre- 
chende Pivotspaltennummer aus einer vorab bereitgestellten Pivotspaltenindizesmatrix ausgelesen wird, 
daB das Nichtvorhandensein einer Pivotspaltennummer zu der voriiegenden Dusenauswahlidentifikations- 
nummer und zu der gewahiten Pivotzeiie durch eInen besonderen Eintrag In der Pivotspaltenindizesmatrix 
eindeutig angezeigt wird, 

und daB die Nichtexistenz einer Losung dadurch angezeigt wird, daB zu kelner moglichen Pivotzeiie in der 
Plvotzeilenindizesmatrix eine Pivotspaltennummer angegeben ist. 

4. Verfahren nach einem der Anspruche 1 bis 3, 

dadurch gekennzeichnet, 
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- daB der Pivotspaltenvektor der Koeffizientenmatrix dss aktuellen Simplextableaus, dessen Werte fOr die Be- 
rechnung neuer Ansteuerwerte fQr die DOsen einer neuen DQsenauswahl benotlgt warden, mit Kenntnis der 
Dusenauswahi vor dam Basisaustauschschritt und der gewahlten Pivotspaltennummer aus elnem vorab be- 
reitgesteilten Feld aller mogiichen optimalen Koeffizientenrnatrizen ausgelesen wird, 

5 - da3 der Vektor der Ansteuerwerte fur die Dusen der Dusenauswalii, die nacli dem Basisaustauschschritt 

vorliegt, berechnet wird und die Werte in diesem Vektor entsprechend der Zuordnung zu den Dusen und 
entsprechend der festgelegten Ordnung der Identifikationsnumnnern fQr die DQsen sortiert werden. 

5. Verfahren nach einem der AnsprQche 1 bis 4, 
10 dadurch gekennzefchnet, 

- da3 mit Kenntnis der Dusenauswahlidentifikationsnummer vor dem Basisaustauschschritt und der gewShlten 
Pivotzeilennummer die Dusenauswahlidentifikationsnummer derjenigen Dusenauswahi, die nach dem Ba- 
sisaustauschschritt voriiegt, aus einer Dusenauswahlidentifikationsnummemmatrix ausgelesen wird 

IS - und da3 nach Vorliegen eines Vektors der Ansteuerwerte, der nur noch nicht negative Werte enthalt, mit 

Kenntnis der dazugehorigen Dusenauswahlidentifikationsnummerder Vektor der Dusenindizes aus einer Ba- 
sisvariablenindizesmatrix ausgelesen wird und die mit ihrer Dusen identifikationsnummer indizierten Dusen 
den entsprechenden Wert aus dem Vektor der Ansteuerwerte (b) erhalten sowie die ubrigen Dusen den An- 
steuenvert Null zugewiesen bekommen und damit der Ansteuen/ektor (k^pt) bestimmt ist. 

20 

6. Verfahren nach elnem der Anspruche 1 bis 5, 

dadurch gekennzeichnet, 

da3 nur diejenigen Spalten von Koeffizientenrnatrizen dual zulassiger Simplextableaus in einer Matrix aller rele- 
vanten Pivotspaitenvektoren abgelegt werden, die bei der Abarbeitung des Verfahrens tatsachlich als Pivotspal- 
25 tenvektoren ausgewahit werden konnen, 

und da3 identische relevante Pivotspaitenvektoren nur einmal in dieser Matrix abgelegt werden. 

7. Vertahren nach Anspruch 6. 

dadurch gekennzelchnet, 

30 

- daS fur die Adressierung der relevanten Pivotspaitenvektoren jedem Pivotspaltenvektor In der Matrix aller 
relevanten Pivotspaitenvektoren eine Pivotspaltenidentlfikatlonsnummer zugewiesen wird und diese in einer 
Pivotspaltenidentifikationsnummernmatrix derart eingetragen werden, daB bel Kenntnis der aktuellen DQsen- 
auswahlidentifikationsnummer und der aktueil gewahlten Pivotzeilennummer die Pivotspaltenidentiflkations- 

35 nummer des aktueil benotigten Pivotspaltenvektors auslesbar ist, 

da3 das Nichtvorhandensein einer PIvotspalte zu einer gewahlten Pivotzeile durch einen besonderen Eintrag 
in der Pivotspaltenidentifikationsnummemmatrix eindeutig angezeigt wird, 

- und da3 die Nichtexistenz einer L6sung dadurch bestimmt wird, daB zu keiner mdglichen Pivotzeile in der 
Pivotspaltenidentifikationsnummemmatrix eine Pivotspaltenidentlfikationsnummer angegeben ist. 

40 

8. Verfahren nach elnem der AnsprQche 3 bis 7, 

dadurch gekennzelchnet, 

daB statt der relevanten Pivotspaitenvektoren nur bestimmte Pivotspaitenvektoren als Reprasentantenvektoren 
im Rechnerspeicher abgelegt werden. wobei die Reprasentantenvektoren eine moglichst kleine Anzahl von Pivot- 

45 spaltenvektoren darstellen, aus denen die tatsachlich benotigten relevanten Pivotspaitenvektoren ermittelt werden, 

und daB jeder Reprasenta.ntenvektor eine Reprasentantenidentifikationsnummer zugewiesen bekommt und zur 
Adressierung der Reprasentantenvektoren eine Reprasentantenidentlfikationsnummernmatrixverwendet wird, in 
derfurjede Diisenauswahlidentlfikationsnummerund jede Pivotzeilennummer die Reprasentantenidentifikations- 
nunnmer desjenlgen ReprSsentantenvektors eingetragen 1st, der fQr die Berechnung des erforderllchen Plvotspal- 

50 tenvektors veiwendet wird. 

9. Verfahren nach Anspruch 8, 

dadurch gekennzeichnet, 

da3 in einer Austauschmatrix fur jede Dusenauswahlidentifikationsnummer und jede Pivotzeilennummer der Zei- 
55 lenindex der auszutauschenden Basisvariablen als Austausch element angegeben ist, um aus dem Reprasentan- 

tenvektor den benotigten Pivotspaltenvektor zu berechnen. 

10. Verfahren nach einem der Anspruche 1 bis 9, 
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dadurch gekennzelchnet, 

daB fQr die Mdglichkeft der Verinderung der DQsenanordnung durch zu erwartende DQsenausf§lle die entspre- 
chenden Daten vorab berechnet und gespetchert warden und nach Feststellen eines Dusenausfails automatlsch 
durch den Rechner auf Verwendung der entsprechenden neuen Daten unngesclialtet wind. 

5 

11. Verfahren nach Anspruchi 10, 

dadurch gelcennzelchnet, 

daB die nach einem Dusenausfall erforderllchen neuen Daten durch eine entsprechende Korrektur und Al<tuali- 
siemng der bis zum DQsenausfail venvendeten Daten bereitgestellt werden. 

10 

12. Verfahren nach einem der AnsprQche 1 bis 11, 

dadurch gekennzeichnet, 

daB das Verfahren unabhangig von der Dimension des Konnmandovektors (bj^om) angewendet wird und dabei fur 
die Berechnung der vorab bereitgesteiiten Daten die nichtauszuiibenden Konnponenten aus der vollstSndigen Sy- 
^5 stemmatrix der DQsenanordnung eiiminlert werden. 



Claims 

20 1 . Method of overall fuel-minimal computer-assisted drive control of a desired determined number and desired de- 
tenmined arrangement of nozzles of a space vehicle, characterised in that 

- a fuel-minimal drive control vector (k^pt) is fonmed with exclusively non-negative drive control values for all 
predetemiined nozzles with the assistance of the dual simplex algorithm, 

25 - for that purpose a nozzle selection identification number (1) is associated with each nozzle selection from the 

total number of nozzles of the space vehicle, which selection represents the optimum nozzle selection for the 
fuel-minimal realisation of a desired command vector (bj^o^^), 

- for ascertaining the optimum drive control, initially the command vector (b^^) containing the force and/or 
moment components to be exerted is transformed Into a start vector (b =bstart) multiplication by a transfor- 

30 matlon matrix ((Bgtart)'^)' 

- thereafter the dual simplex algorithm Is performed until the vector (b) of the drive control values contains only 
non-negative drive control values or, however, the non-existence of a solution is Indicated, 

- and In that case the dual simplex algorithm is realised In that firstly all data necessary for the processing of 
the dual simplex algorithm and able to be computed off-line beforehand are ascertained, are filed in data fields 

35 In the computer memory and on use of the dual simplex algorithm during on-line operation are called up from 

the computer store by appropriate addressing and employed in the computation, wherein processing is un- 
dertaken only with a reduced simplex tableau in that merely the nozzle selection identification number and the 
associated vector (b) of the drive control values are recorded. 

40 2. Method according to claim 1, characterised in that the transfomiatlon matrix ((Bstart)'^) represents the inverted 
basic matrix of those nozzles which belong to the starting nozzle selection, this nozzle selection is optimal for a 
desired command vector, a starting nozzle selection identification number (i^tart) associated with this nozzle 
selection and the start vector (b ^ b^tait) contains the drive control values of these nozzles. 

<5 3. Method according to one of claims 1 and 2, characterised in that 

each nozzle of the space vehicle receives an identification number and the values in the vector of the drive 
control values for a specific nozzle selection are always classified in a fixed order on the basis of the identifi- 
cation numbers of the nozzles, 

so . in the case of processing the dual simplex algorithm after the ascertaining of a pivot line number and with 
knowledge of the nozzle selection Identifteation number of the nozzle selection actually present the corre- 
sponding pivot column number is read out of a pivot column indice matrix prepared beforehand, 
the non-presence of a pivot column number with respect to the nozzle selection identification number present 
and to the selected pivot line is unambiguously indicated by a special entry in the pivot column indice matrix 

55 . and the non-existence of a solution is indicated in the manner that a pivot column-fHimberHSHndicated-with 

respect to no possible pWot line in the pivot line Indice matrix. 

4. Method according to one of claims 1 to 3, characterised in that 
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- the pivot column vector of the coefficient matrix of the actual simplex tableau of which the values are required 
for the computation of new drive control values for the nozzles of a new nozzle selection is read out of a 
previously prepared field of alt possible optimum coefficient matrices with knowledge of the nozzle selection 
before the basic exchange step and of the selected pivot column number 

5 - and the vector of the drive control values forthe nozzles of the nozzle selection present afterthe basic exchange 

step Is computed and the values in this vector are classified in correspondence with the association with the 
nozzles and in correspondence with the fixed association of the identification numbers forthe nozzles. 

5. Method according to one of claims 1 to 4, characterised in that 

10 

- the nozzle section identification number of that nozzle section which is present after the basic exchange step 
is read out of a nozzle selection identification number matrix with Icnowledge of the nozzle selection identifi- 
cation number before the basic exchange step and of the selected pivot line number 

and after presence of a vector of the drive control values which contain only non -negative values the vector 
15 of the nozzle indice is read out of a basic variable in dice matrix with l<n owl edge of the associated nozzles 

selection identification number, and the nozzles denoted by their nozzle identification numbers receive the 
corresponding value from the vector (b) of the drive control values and the remaining nozzles are assigned 
the drive control value zero and thus the drive control vector (j^pt) is detennined. 

20 6. Method according to one of claims 1 to 5, characterised in that only those columns of coefficient matrices of a dual 
permissible simplex tableau are filed in a matrix of all relevant column vectors which can be actually selected as 
pivot column vectors in the processing of the method and that identical relevant pivot column vectors are filed only 
once in this matrix. 

25 7, Method according to claim 6, characterised in that 

for the addressing of the relevant pivot column vectors a pivot column identification number is assigned to 
each pivot column vector in the matrix of all relevant pivot column vectors and these are entered in a pivot 
column identification number matrix in such a manner that the pivot column identification number of the actually 
30 required pivot column vector can be read out with knowledge of the actual nozzle section identification number 

and the actually selected pivot line number, 

the non-presence of a pivot column with respect to a selected pivot line is unambiguously indicated by a special 
entry in the pivot column identification number matrix 

and the non-existence of a solution is detemiined in the manner that a pivot column identification number is 
35 indicated with respect to no possible, pivot line in the pivot column identification number matrix. 

8. Method according to one of claims 3 to 7, characterised in that instead of the relevant pivot column vectors only 
specific column vectors are filed in the computer memory as representative vectors, wherein the representative 
vectors represent a smallest possible number of pivot column vectors from which the actually required relevant 

40 pivot column vectors are ascertained and that each representative vector is assigned a representative identification 

number and a representative identification number matrix is used forthe addressing of the representative vectors, 
in which matrix for each nozzle selection identification number and each pivot line number the representative 
identification number of that representative vector is entered which is used for the computation of the required 
pivot column vector. 

45 

9. Method according to claim 8. characterised in that the line index of the basic variables to be exchanged is Indicated 
as an exchange element in an exchange matrix for each nozzle selection identification number and each ph^ot line 
number in order to compute the required pivot column vector from the representative vector. 

50 1 0. Method according to one of claims 1 to 9, characterised in that for the possibility of variation in the nozzle anBnge- 
ment by anticipated nozzle failures the con^espondlng data are computed beforehand and stored and after deter- 
mining a nozzle failure there is automatic switching over by the computer to use of the appropriate new data. 

11. Method according to claim 1 0, characterised in that the new data required after a nozzle failure are prepared by 
55 an appropriate correction and actualisation of the data used up to the nozzle failure. 

12. Method according to one of claims 1 to 11 , characterised that the method is used independently of the dimension 
of the command vector f^m) ^^^^ c^®® computation of the previously prepared data the components 
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which are not to be exerted are eliminated from the complete system matrix of the nozzle arrangement 
Revendications 

5 

1. Proc6d6 pour ia commande globale asslst§e par calculateur, corresponolant k une consommatlon minimale de 
propergol, d'un nombre et d'une disposition fix^ & volenti de tuyeres sur un v^hlcule spatial, 
caracteris^ en ce 

10 - qu'un vecteur de commande (kgpt) con^espondant k une quantity minimale de propergol et poss^dant des 
valeurs de commande exclusivement non negatives pour toutes les tuyeres pr^determindes est fonmS k I'alde 
de Talgorithme du simplexe dual, 

qu'a cet effet un numero (I) d'identification de selection de tuydres est associe k chaque selection de tuyeres 
pamii latotalite de tuyeres du vehlcule spatial qui represente la selection optimale de tuyeres pour la realisa- 
15 tion, correspondant k une consommatlon minimale de propergol, d'un vecteur quelconque de commande 

que pour la d6termination de la commande optimale, le vecteur de commande i^n), qui contient la compo- 
sante de force devant §tre appliqute et/ou la composante de moment devant etre appilquee. est transfomn^ 
par multiplication par une matricede transformation ((Bd^manage)'^) ®" un vecteur de ddman'age (b=bd6marrage). 
20 - qu'ensuite I'algorithme du simplexe dual est exteutS Jusqu'^ ce que le vecteur d'une valeur de commande (b) 
contienne unlquement des valeurs de commande non negatives ou bien que la non existence d'une solution 
soit indiqu^e, et 

que Talgorithme du simplexe dual est mis en oeuvre par le fait qu'au prealable toutes les donnees, qui sont 
necessaires pour letraitement de I'algorithme du simplexe dual peuvent etre calculees off-line (au prealable), 
25 sont d6tennin6es, sont m6moris§es dans des champs de donn6es dans la memoire du calculateur et sont 

appelees k partlr de la memoire du calculateur lors de I'utilisation de Talgorithme du simplexe dual pendant 
le fonctionnement on-line au moyen d'adressages appropri^s et sont utilises dans le calcul, le travail etant 
effectud uniquement avec un tableau slmpllcial r^dult, par le fait que seul le numero d'identification de selection 
de tuy&res et le vecteur, qui est associ^, des valeurs' de commande (b) sont pilot^s. 

30 ... 

2. Proc6d6 seion la revendication 1 , caract§ris6 en ce 

que la matrice de transformation ({Bdemarrage)"^) represente la matrice de base inversee destuyferes, qui font partie 
de la selection des tuyeres de d6marrage, que cette selection de tuyeres est optimale pour un quelconque vecteur 
de commande, qu'un numero (idemarrage) d'identification de selection de tuyeres est associe k cette selection de 
35 tuyeres et que le vecteur de d6marrage (bsbd^manage) contient les valeurs de commande de ces tuyeres. 

3. ProcSdd selon la revendication 1 ou 2, caracteris^ en ce 

- que chaque tuy6re du v6hicule spatial regoft un num§ro tfidentiflcation, et que les valeurs contenues dans le 
40 vecteur des valeurs de commande pour une s61ection d6tenDin6e des tuyferes sont en pemrianence classees 

selon un ordre fix6, en r6f6rence aux num6ros d'identification des tuyferes, " 

- que lors de I'ex^cutlon de I'algorithme du simplexe dual, apres la d6terminatlon du num6ro de llgne-pivot avec 
la connalssance du num6ro d'identification de la selection actuelle des tuyferes, le num6ro correspondant de 
colonne-pivot est lu a partir d'une matrice d'indices de colonnes pivots foumies, 

45 . que I'absence d'un numero de colonne-pivot en rapport avec le numero present d'identification de selection 

de tuyeres et avec la llgne-pivot selection n6e soit affich6e de fafon nette au moyen d'une entree particuliere 
dans la matrice d'indices de colonnes-pivots, et 

- que la non existence d'une solution est indiqu6e par le fait qu'aucun num6ro de colonne-pivot n'est Indiqu6 
pour aucune llgne-pivot possible dans la matrice d'indices de llgnes-pivots. 

50 

4. ProcSdS selon I'une des revendications ^ k3, caract^rls^ en ce 

- que le vecteur de colonne-pivot de la matrice de coefficients du tableau simpllcial actual, dont les valeurs sont 
necessaires pour le calcul de nouvelles valeurs de commande pour les tuyeres d'une nouvelle selection de 

55 tuyeres, est lu, sur la base de la connaissance de la selection des tuyeres avant I'^tape d'6change de base 

et avec la connalssance du num§ro de colonne-pivot s^lectlonnd k partir d'un champ, pr6par§ au prSalable, 
de toutes les matrices possibles de coefficients, 

- que le vecteur des valeurs de commande pour les tuyeres de la selection de tuyeres, qui est present apres 



15 



V i 



EP 0 977 687 B1 

I'dtape d'^hange de base, est calculi et que les valeurs contenues dans ce vecteur sont class^es en fonctlon 
de rassociation aux tuy&res et conformSment d I'ordre fixS des numSros d'identification pour les tuyeres. 

5. Proc6d§ selon I'une des revendications 1 k 4, caract6ris6 en ce 

5 

qu'avec la connalssance du numero d'identification de selection de tuyeres avant I'etape d'^change de base 
et avec la connalssance du numero s§iectlonn6 de llgne-plvot. le numero d'identification de la selection de 
tuyeres, qui est present apr&s I'etape d'^hange de base, est lu dans une matrice de num^ros d'identification 
de selection de tuy&res, et 

10 - qu'apr^s I'apparitlon d'un vecteur des valeurs de commande. qui contient encore seulement des valeurs non 
negatives, et avec la connalssance des numeros d'identification de selection de tuyeres, qui y sont associes, 
le vecteur des indices des tuyeres est lu k partir d'une matrice d'indices de variables de base et les tuyeres, 
qui sont indicees avec leurs numeros d'identification de tuydres, regoivent la valeurcorrespondante k partir 
du vecteur des valeurs de commande (b), et que la valeur de commande z6ro est affectde aux autres tuyferes 

IS et que par consequent le vecteur de commande i^^) est d6termin6. 

6. Precede selon I'une des revendications 1 k 5, caractSris^ en ce . .. 

que seules les colonnes de matrices de coefficients de tableaux simpliciaux admissibles sont memorisees 
20 dans une matrice de tous les vecteurs importants de coionnes-pivots, qui peuvent Stre effectivement sSiec- 

tionnes en tant que vecteurs de coionnes-pivots, lors de la mise en oeuvre du proc6d6, et 
que des vecteurs identiques importants de coionnes-pivots sont m^morisdes seulement une fols dans cette 
matrice. 

25 7. Proc§d6 selon la revendication 6, caractSrisS en ce 

que pour I'adressage des vecteurs importants de coionnes-pivots, un numero d'identification de colonne-pivot 
est affects k chaque vecteur de colonne-pivot de la matrice de tous les vecteurs importants de coionnes-pivots 
et que ce numSro est Introduit dans une matrice de numSros d'identification de coionnes-pivots de telle sorte 
30 que, moyennant la connaissance du num6ro actual d'identification de selection de colonne et du num6ro 

actueilement s6lectionn§ de ligne-pivot, le numero d'identification de colonne-pivot du vecteur actuellement 
n6cessaire de coionnes-pivots peut dtre lu, 

que I'absence d'une colonne-pivot en rapport avec une ligne-pivot s6lectionn6 est affich6 de fafon nette au 
moyen d'une entr6e particuliere dans la matrice de numeros d'identification d'une colonne-plvot, et 
35 - que I'absence d'une solution est dStemiinSe par le fait qu'aucun numSro d'identification de colonne-pivot n'est 

indiquS pour aucune ligne possible-pivot dans la matrice de numeros d'identification de coionnes-pivots. 

8. ProcSdS selon Tune des revendications 3 d 7, caractSrisg en ce 

40 qu'k la lace des vecteurs importants de coionnes-pivots, seuls certains vecteurs de coionnes-pivots sont me- 

morises en tant que vecteurs repr§sentatlfs dans la m^moire du calculateur, les vecteurs representatifs re- 
pr^sentant un nombre aussi faible que possible de vecteurs de coionnes-pivots, k partir desquels sont deter- 
mines les vecteurs importants de coionnes-pivots, qui sont effectivement n6cessaire, et 
qu'un numero d'identification de reprdsentant est affecte a chaque vecteur representatif et que pour I'adres- 

45 sage des vecteurs representatifs on utilise une matrice de numeros d'identification de representants, dans 

laqueile est produrt, pour chaque numero d'identification de selection de tuyere et pour chaque numero de 
llgne-prvot, le numero d'identification de representant du vecteur representatif. qui est utilise pour le calcul du 
vecteur necessaire de colonne-pivot. 

50 9. Procede selon ia reyendication 8, caracterise en ce que 

dans la matrice d'echange, pour chaque numero d'identification de selection de tuyeres et pour chaque hum6ro 
de ligne-pivot, I'indice de ligne de base devant etre echange est indique en tant qu'eiement d'echange, pour le 
calcul du vecteur necessaire de coionnes-pivots k partir du vecteur representant. 

55 10. Precede selon I'une des revendications 1 k 9. caracterise en ce 

que pour la possibillte de modifier la disposition des tuyeres en ralson de pannes de tuyeres auxquelles 11 faut 
s'attendre, les donnees correspondantes sont calcuiees au prSalable et memorisees et, apres la determination 
d'une panne de tuyere, une commutation automatique est effectuee par le calculatebr, sur rutilisatlon des nouvelles 
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donn6es correspondantes. 

1 1 . Proc^de selon ia revendication 1 0, caracterisd en ce 

que les nouvelles donn^es n^essalres apr^s une panne de tuyere sont pr^ardes an moyen d'une correction et 
d'une actualisation correspondantes des donnSes utilisabies jusqu'^ la panne d'une tuyere. 

12. Proc6d6 selon Tune des revendications 1 & 11, caract6ris§ en ce 

que le proc6d6 est utilise ind^pendamment de la dimension du vecteur de commande (b^n,) et que. pour le calcul 
des donn^es pr6par6es au pr6alable, les composantes, qui ne doivent pas §tre actlv6es, sont 6limin6es de la 
matrice complete du syst&me de I'ensemble des tuyeres. 
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